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In this paper we formulate the theory of the interaction of a diatomic hnear molecule in a spatially 
degenerate state with the non-resonant laser field and of the rovibrational dynamics in the presence 
of the field. We report on ab initio calculations employing the double electron attachment inter- 
mediate Hamiltonian Fock space coupled cluster method restricted to single and double excitations 
for all electronic states of the Rb2 molecule up to 5s -I- 5d dissociation limit of about 26.000 cm~^. 
In order to correctly predict the spectroscopic behavior of Rb2, we have also calculated the electric 
transition dipole moments, non-adiabatic coupling and spin-orbit coupling matrix elements, and 
static dipole polarizabilities, using the multireference configuration interaction method. When a 
molecule is exposed to a strong non-resonant light, its rovibrational levels get hybridized. We study 
the spectroscopic signatures of this effect for transitions between the X^E^ electronic ground state 
and the A^Ej and b^n„ excited state manifold. The latter is characterized by strong perturba- 
tions due to the spin-orbit interaction. We find that for non-resonant field strengths of the order 
lO^W/cm^, the spin-orbit interaction and coupling to the non-resonant field become comparable. 
The non-resonant field can then be used to control the singlet-triplet character of a rovibrational 
level. 



I. INTRODUCTION 

Rubidium was one of the first species to be Bose- 
condensed and nowadays it can routinely be cooled 
and trapped. It has therefore become the drosophila 
of ultracold physics. Its long-range interatomic interac- 
tions have extensively been studied, and this has allowed 
to very accurately determine the scattering length and 
Cg coefficient @t3|- ^^^2 molecules have been formed 
out of ultracold rubidium atoms using both photo- and 
magneto-association 0]. Photoassociation and Fes- 
hbach spectroscopy have also served to measure the 
low- lying shape resonances of the rubidium dimer 0411 ■ 
Trapping rubidium in an optical lattice has facilitated 
studies of atom-molecule dark states and transfer- 
ring the molecules into their vibrational ground state [ll| . 
The Rb2 molecule continues to draw attention in the con- 
text of the coherent control of ultracold collisions [T2| - [lH 
and femtosecond photoassociation ^16|-[T9j. These exper- 
iments as well as those e mployi ng photoassociation with 
continuous wave lasers |2Cll424| require precise spectro- 
scopic knowledge not only of the ground but also the 
excited states for both interpretation and detection. 

The electronic ground and excited states have ex- 
tensively been studied. According to Huber and 
Herzberg [2^, the Rb2 molecule was first observed in 
a spectroscopic experiment by Lawrence and Edlefsen as 
early as 1929 [1^. Since then, dozens of experimental 
papers have appeared, and recently the Rb2 molecule 
has been revisited in the context of cold molecules stud- 
ies. The ground X^S+ state has been investigated in 



Ref. [27|, while the most accurate experimental results 
for the a^S+ state have been reported by Lozeille et 
al. [28l. Beser et al. [11] and Tiemann and collabora- 
tors [30| . The most important excited states correspond- 
ing to the ^S-I-^P dissociation limit, the A^E+ and b^n„ 
states, have extensively been analysed in Ref. Less 
experimental information is available for other excited 
states. Notably, the (l)'^Sj" state has been studied in 
Ref. [sij, and Ref. (s^j reports the experimental obser- 
vation of the (2)^ng state. The pure long-range state of 
0^ symmetry, that is important for the photoassociation 
of ultracold Rb atoms, has been analysed in Ref. (33| . 
Several of these experimental data were successfully em- 
ployed to derive empirical potentials that reproduce the 
spectroscorac data with the experimental accuracy, cf. 
Refs. [13,1311 for the ground state X^T.+ and Refs. [11, 
for the a'^E^ potential. The coupled manifold of the 
A^E+ and b^n^ states was deperturbed by Bergeman 
and collaborators [H, [s^ with the corresponding po- 
tential energy curves and spin-orbit coupling matrix ele- 
ments reported in Ref. [33 |. Potential energy curves for 
other electronic states fitted to the experimental data 
are older, cf. Ref. [Ill for the empirical potential energy 
curve of the (l)"^ng state, and Refs. [11] and [13] for those 
of the (2)^E+ state and {2Yllu states, respectively. 

Given this extensive amount of experimental data, it 
is not surprising that many theoretical calculations have 
tackled the ground and excited states of the rubidium 
dimer. The first ab initio calculation on the Rb2 molecule 
dates back to 1980 and was reported by Konowalow and 
Rosenkrantz [11]. In the 1980's and 1990's several elec- 
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tronic structure calculations of different quality followed. 
Three more recent studies reported ab initio data of vary- 
ing accuracy for the potential energy curves and in some 
cases other properties such as couplings and transition 
moments of Rb2. In 2001, Park et al. [39J reported 
non-relativistic potentials for all molecular states. The 
root mean square deviation (RMSD) between the theo- 
retical well depths and the available experimental data 
was 235 cm~^, i.e., 9.9% on the average. The 2003 cal- 
culations by Edvardsson et al. [i^l were devoted to the 
ground state potential and six excited state potentials of 
ungerade symmetry. The spin-orbit coupling matrix ele- 
ments were also reported. The overall accuracy of these 
results was about the same as in Ref. [s^ with a RMSD 
of 180 cm~^ representing an average error of 25%. Note 
that since the number of states considered in Refs. [1^ 
and [i^l differs, the absolute RMSD may be smaller and 
the percentage error larger. Finally, in 2012 AUouche and 
Aubert-Frecon ^] reported calculations of all molecu- 
lar states and spin-orbit coupling matrix elements corre- 
sponding to the dissociation limits 5s -|- 5s, 5s -I- bp, and 
5s -|- Ad. These calculations are much more accurate than 
any other previously reported in the literature with a 
RMSD of 129 cm~^, i.e., an error of 5.5% only. However, 
they do not cover highly excited molecular states that are 
of interest for the detection of ultracold molecules [42| as 
well as photoassociation into states with ion-pair charac- 
ter [ist - fisj . 

Photoassociation into highly excited electronic states 
is at the core of a recent proposal for the production of ul- 
tracold Rb2 molecules [45| , aimed at improving the initial 
femtosecond experiments [T6l - [l9| . It is based on multi- 
photon transitions that can easily be driven by femtosec- 
ond laser pulses and allow to fully take advantage of the 
broad bandwidth of femtosecond laser pulses while driv- 
ing the narrow photoassociation transition j46i] . More- 
over, multi-photon photoassociation populates highly ex- 
cited electronic states with ion-pair character and strong 
spin-orbit interaction. This features are advantageous for 
an efficient stabilization of the photoassociated molecules 
into deeply bound molecules in the electronic ground 
state 4^ . The theoretical modeling of the proposed pho- 
toassociation scheme required the knowledge of precise ab 
initio potential energy curves, spin-orbit and nonadia- 
batic coupling matrix elements, electric transition dipole 
moments and dynamic Stark shifts. All these data were 
not available in the literature. Moreover, the newly de- 
veloped tools of electronic structure theory based on the 
Fock space coupled cluster method [47'J49'] could possibly 
allow for reaching a better accuracy of the potentials than 
reported in Refs. [39l - l4l| . Last but not least, calculations 
of the electric properties for diatomic molecules in spa- 
tially degenerate electronic states are scarce. To the best 
of our knowledge, only two studies considered this prob- 
lem, in the context of the dispersion interactions between 
molecules (sol . [sij rather than non-resonant interactions 
with an external field, and a systematic theoretical ap- 
proach has not yet been proposed. Moreover, the pres- 



ence of spin-orbit coupling between the electronic states 
has been neglected in a recent treatment of nuclear dy- 
namics in a non-resonant field [s^ [ssj . Such an approxi- 
mation does not allow to study the competition between 
the spin-orbit coupling and the interaction with a non- 
resonant field which may both significantly perturb the 
spectrum. 



Here, we fill this gap and report a theoretical formu- 
lation of the interaction of a H state molecule in a non- 
resonant field and of the rovibrational dynamics in the 
presence of the field. We also report ab initio calcula- 
tions of all potential energy curves, spin-orbit and nona- 
diabatic coupling matrix elements corresponding to the 
dissociation limits up to and including 5s -I- 5d. We test 
our ab initio results by comparing the main spectroscopic 
characteristics of the potentials to the available experi- 
mental data. We devote a special emphasis to the impor- 
tant manifold of the A^S+ and b^H^ states, comparing 
our results to Refs. [2O, [s^. Since the electric properties 
of spatially degenerate electronic states were not exten- 
sively studied in the literature thus far, we report here 
first, to the best of our knowledge, ab initio calculations 
of the irreducible components of the polarizability tensor, 
including their dependence on the interatomic distance 
R, for the A^E+ and b'^Hu states. Finally, we study the 
effect of a non-resonant field on the spectroscopy in the 
A^I]+ and b'^H„ manifold. This is motivated by our re- 
cent proposal for enhancing photoassociation by control- 
ling shape resonances with non-resonant light js^ . [ssj . In 
order to significantly modify the scattering continuum of 
the atom pairs to be photoassociated, rather large non- 
resonant intensities are required. Since the bound rovi- 
brational levels are much more affected by a strong non- 
resonant field than continuum states, it is important to 
investigate how the corresponding spectroscopic features 
change. 



Our paper is organized as follows. In Sec. |TT] we for- 
mulate the theory of the interaction of a homonuclear 
molecule with an external non-resonant field. In Sec. IIIII 
we provide the theoretical description of the perturbation 
of spectra by a non-resonant field, using as an example 
the spin-orbit coupled manifold of the A^I]+ and b^n„ 
electronic states of Rb2 . We briefly summarize the ab ini- 
tio methods employed in our calculations in Sec. IIVI and 
discuss the results of these calculations in Sec. [V] In par- 
ticular, we compare our data with results available in the 
literature and discuss the ability of the ab initio results 
to reproduce the high-resolution spectroscopic data for 
the A^E+ and b^H„ manifold [20,|3i|. We then describe 
the interaction with a non-resonant field and study its 
spectroscopy signatures on the transitions between the 
electronic ground state and the A^I]+ and b^H,( mani- 
fold. Finally, Sec. IVII concludes our paper. 



3 



II. LINEAR MOLECULE IN A 
NON-RESONANT ELECTRIC FIELD 

We consider the interaction of a linear diatomic 
molecule with an electric field with the direction taken 
along the Z axis of the space-fixed coordinate system, 
E = (OjO,!?). To the second order, the Hamiltonian for 
the interaction of the molecule with the electric field can 
be written as. 
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int 



(1) 



where d^^ and denote the appropriate components 
of the electric dipole moment and electric dipole polar- 
izability in the space-fixed frame. Since we deal with 
a homonuclear molecule, only the second term of the 
above Hamiltionian will be relevant in the present analy- 
sis. To evaluate the matrix elements of the Hamiltonian 
in the electronic and rovibrational basis, we rewrite a^^^ 
in terms of the polarizability components in the body- 
fixed frame. The dipole polarizability component 
can be expressed in terms of space-fixed irreducible ten- 
sor components Ofm'^^ 
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For the irreducible tensor components, the transforma- 
tion from the space-fixed to the body-fixed coordinate 
system is given by the rotation matrices f,{R), 



(i),BF 
k 



(3) 



k=-l 



Hence, we have 
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For simplicity, we omit the superscripts SF/BF in the 
rest of the paper as from now we will use only the body- 
fixed quantities. For a diatomic linear molecule the set 
of the Euler angles R can be chosen as i? = (0,0,0), 
where 6 is the angle between the molecular axis and the 
space-fixed Z axis. We assume in this paper that the 
molecular axis defines the body-fixed z axis. Note also 
that for such a specific choice of the Euler angles, the 
Wigner D functions appearing in Eqs. (|4]) reduce to: 



{l-k)\ 
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(5) 



where Pj^ are the associated Legendre polynomials. For 
any linear molecule, the non-zero irreducible components 



of the dipole polarizability are o'q^ and a)^' ■ In addition, 
for a linear molecule in a H electronic state, the and 
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(2) 

terms do not vanish. They should be viewed as 
off-diagonal polarizability tensor components connecting 
two degenerate electronic states, |Hi) and |H_i), with 
opposite projection of the total electronic orbital angular 
momentum on the molecular axis. See, for instance, Eq. 
(16) of Ref. 151]. 

The non-vanishing body-fixed polarizability compo- 
nents are most conveniently expressed in terms of the 



Cartesian tensor elements aa^ i = x,y,z. 
related to the trace of the polarizability. 



Then al^"^ is 
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Q!q to the anisotropy of the polarizability. 
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and, for a molecule in a H electronic state, a_2 and 
reflect the difference between two perpendicular compo- 
nents, 



(2) (2) 

Q^2 = ^-2 = %y ^ ^xx- 



(8) 



For a linear molecule in a S state, the definitions of 
the Cartesian components of the polarizability tensor an 
in Eqs. ([6]) to (l8|) are unambigous. The zz and xx com- 
ponents are simply the parallel and perpendicular com- 
ponents, a|| and a±, respectively. Thus, the irreducible 
tensor components appearing in Eqs. (I6|) to ([8]) are sim- 
ply related to the trace a and the anisotropy Aa of the 
polarizability tensor, 



r = ~V3a, 
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(9) 



Obviously, for a S state molecule the xx and yy compo- 
nents are equal, and aj^'' — 0. 

In the case of a molecule in a degenerate electronic 
state (n, A etc.) some caution is needed when employ- 
ing the Cartesian components an, since one has to specify 
the basis of the electronic states, in which these quan- 
tities are expressed. Equation ^ assumes the Carte- 
sian components, ayy and a^x, to be calculated for the 
|Hj;) state. However, the Cartesian basis {|Ha;), |Hj,)} for 
the H electronic state is not convenient for the dynami- 
cal calculations, since the spin-orbit coupling matrix el- 
ements are complex in this basis. Therefore, we prefer 
to use the spherical basis {|H_i), |Hi)} for the H state 
over the Cartesian basis {jH^;), |Hy)} since it avoids com- 
plex quantities in the calculations and allows for a sim- 
ple adaptation of the Hund's case (a) wave function to 
a given symmetry of the rovibrational level. Therefore, 
we will use the irreducible polarizability components am 
rather than the Cartesian an. 

Combining Eqs. ([T|) to (|S]) and making use of proper- 
ties of the rotation matrices f,{R), one arrives at the 
following Hamiltonian for the interaction of the homonu- 
clear diatomic molecule with the static electric field, 



4 



H\„t = — — 



£2 

2 



T3 



(2) 



P2°(C0S( 



, 1 (2) 



P2 (cos( 



+ 4aL1P2" 



'(cos 61) 



(10) 



r 



The above Hamiltonian is valid for any isolated elec- 
tronic state of a linear homonuclear molecule. Albeit, 
the last two terms in this equation are relevant only for 
molecules in a 11 electronic state. Let us stress here that 
although this form of the Hamiltonian seems a bit elab- 
orate at first glance, it simplifies the evaluation of the 
matrix elements in the symmetry-adapted basis set, and 
it also avoids any ambiguities when employing the Carte- 
sian polarizability components for degenerate electronic 
states. Equation (fTO|) also assumes the frequency of the 
non-resonant field to be far from any resonance and the 
two-photon rotating-wave approximation. 



III. HAMILTONIAN FOR THE RB2 MOLECULE 
IN THE MANIFOLD OF THE COUPLED A^E+ 
AND B=*n„ EXCITED STATES INTERACTING 
WITH A NON-RESONANT FIELD 



basis functions |n, A)\S, S) | J, il, M) that are products of 
the electronic |n,A), electron spin [S*, S) and rotational 
\J,n,AI) functions. Here, J is the total angular mo- 
mentum quantum number, S is the total electronic spin 
quantum number, A and E are the projections of the total 
electronic orbital and total electronic spin angular mo- 
menta onto the molecular axis, and M is the projection 
of the total angular momentum onto the Z space-fixed 
axis, n labels the nonrelativistic dissociation limit of the 
molecular state. We also define the projection of the to- 
tal, electronic orbital plus spin, angular momentum onto 
the molecular axis, i7 = A -I- S. For the coupled A^I]+ 
and h'^Hu manifold, we consider the rovibrational lev- 
els of the e spectroscopic symmetry and odd parity. For 
simplicity, any hyperfine structure effects are neglected 
here. The properly symmetry-adapted Hund's case (a) 
wavefunctions read, 



We construct the Hamiltonian for the nuclear motion 
in Hund's case (a) coupling scheme with the primitive 



J 



|AiSo+,J,M,e) = |A,0)|0,0)|J,0,Af) , 
|b3Ho+,J,M,e 

|b3H2„, J,M, e 



^[|b,l)|l,-l)|J,0,M)-|b,-l)|l,l)|J,0,M)] , 



(11) 



V2 



[|b,l)|l,l)|J,2,M)-|b,-l)|l,-l)|J,-2,Af)] 



r 



The first two states have a projection of the total angu- 
lar momentum onto the molecular axis \Q\ — 0, while the 
third one has jfij = 2. In the field- free case, the state with 
\n\ = 2 is decoupled from the states with \^\ =0, and it 
is not accessible from the ground electronic state in the 



J 



one-photon dipolar transitions considered here. Conse- 
quently, the field-free model Hamiltonian Hq describing 
the nuclear motion in the manifol of the coupled A^S^ 
and h^Hu states can be represented by following 2x2 
matrix. 



dHn = 



2iiw + y 



(P) 



_j 



r 



(12) 



where T 



and rotational kinetic energy operators and 



is the sum of the vibrational 
V\R), 



k = A^S+j b'^Hu, denotes the respective potential 



energy curves in the Born-Oppenheimer approxima- 
tiom ii{R) = (Ail]+|i?so|b3n„) 101=0 and 6(i?) = 
(b'^H„|i/so|b'^H„)|Q|^o are the spin-orbit coupling ma- 
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trix elements, and only the electronic states with = 
are included. 

When the electric field is switched on, the two compo- 
nents of the h'^Hu state with opposite A projections are 
coupled. Therefore, not only the interaction Hint from 
Eq. ([T0| has to be added to the Hamiltonian Hq for the 



A-'^E+ and b'^n„ states, but also the matrix ([T2|) has to 
be extended so as to include the jfij = 2 component orig- 
inating from the b'^n„ state. Hence, in the presence of 
the electric field the rovibrational levels of the A^I]+ and 
b'^n„ manifold are obtained by diagonalizing the Hamil- 
tonian represented by the following 3x3 matrix. 



J 



H = 







T 



Wo/2{R,0) T 





Wo/2iR,0) 



UR)) 



(13) 



r 



The diagonal elements of the interaction potentials in- 
corporating the interaction with non-resonant field are 
given by. 



W''{R,e) = V^{R)+H^^ 



(14) 



where k =A-^E+ or b'^H„ and ifj^^ is given by Eq. (fTO|) 
for the electronic state labeled by k. The off-diagonal 
term due to the non-resonant field, Wo/2(i?,6'), couples 
the = 0+ and — 2„ components resulting from 
the b'^Hii state. It is proportional to the off-diagonal 
polarizability of the molecule in the b'^H state. 



w^/2{R.e)^-^-K 



(i?)F|(cos0), (15) 



with a\ ' defined by Eq. Analogously to Eqs. ^ 
and (HH), the Hamiltonian for the molecule in its elec- 
tronic ground state interacting with a non-resonant field 
is simply given by T + W^^'^s (i?, 9). 



IV. AB INITIO ELECTRONIC STRUCTURE 
AND DYNAMICAL CALCULATIONS 

We adopt the computational scheme successfully ap- 
plied to the ground and excited states of the calcium 
dimer j[55-59|, magnesium dimer [60l l6ll . strontium 
dimer [62l . |63|, (BaRb)"'' molecular ion [64)|, and SrYb 
heteronuclear molecule ;65i] . The potential energy curves 
for the singlet and triplet gerade and ungerade states of 
the Rb2 molecule corresponding to the first seven lowest 
dissociation limits, 5s -I- 5s, 5s -I- 5p, 5s -I- 4(i, 5s -I- 6s, 
5s -I- 6p, hp + bp, and 5s -I- 5(i, have been obtained by a 
supermolecule method. 



V 



(16) 



where denotes the energy of the dimer computed 

using the supermolecule method (SM), and -E'x'^, X — A 
or B, is the energy of the atom X. The calculations for 
the excited states employed the recently introduced Dou- 
ble Electron Attachment Intermediate Hamiltonian Fock 



Space Coupled Cluster method restricted to single and 
double excitations (DEA-IH-FS-CCSD) [13-111 • Start- 
ing with the closed-shell reference state for the doubly 
ionized molecule Rb2~'' that shows the correct dissocia- 
tion at large interatomic separations, R, into closed-shell 
subsystems, Rb++Rb+, and using the double electron 
attachment operators in the Fock space coupled cluster 
ansatz makes our method size-consistent at any inter- 
atomic separation R and guarantee the correct large- i? 
asymtptotics. Thus, the DEA-IH-FS-CCSD approach 
overcomes the problem of the standard coupled clus- 
ter method restricted to single and double excitations 
(CCSD) and the equation of motion CCSD method [H] 
with the proper dissociation into open-shell atoms. The 
potential energy curves obtained from the ab initio cal- 
culations were smoothly connected at intermediate in- 
teratomic separations with the asymptotic multipole ex- 
pansion [H^l • The Cg coefficient of the electronic ground 
state and the C3 coefficient of the first excited state 
were fixed at their empirical values derived from high- 
resolution spectroscopic experiments [1, [33|. while the 
remaining coefficients were taken from Ref. [66{ . 

The transitions from the ground X^E+ state to the 
and ^H„ states and from the a^E+ to the ^Ej" and 



^ to the ^E+ 

^lig states are electric dipole allowed. The transition 
dipole moments for the electric transitions were com- 
puted from the following expression (67| . 



iv+ 
9 



d^{n ^ a) = ( a^E+ 



{nf\K[ 



(17) 



where r.^, i = x,y or z, denotes the ith component of the 
position vector. Note that in the first term of Eq. p7)) 
i ~ X or y corresponds to transitions to ^H^ states, while 
i — z corresponds to transitions to ^E+ states. The tran- 
sitions from the a^E+ state connect this state with the 
'^Hg and '^E^ states, through the x and y and z operators, 
respectively. 

We expect the rovibrational energy levels of the excited 
electronic states of Rb2 to show perturbations due to the 
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nonadiabatic couplings between the states. Analysing 
the pattern of the potential energy curves, we have found 
that many potential energy curves display avoided cross- 
ings, suggesting strong radial couplings between these 
electronic states. We have therefore computed the most 
important radial coupling matrix elements, defined by 
the expression. 



\g/u 



OR 



(18) 

where n ■(r^ n' signifies that the electronic states n and 
n' are coupled. Note that the radial derivative operator 
couples states with the same projection of the electronic 
orbital angular momentum on the molecular axis A. 

Electric transition dipole moments, radial non- 
adiabatic coupling and spin-orbit coupling matrix ele- 
ments were obtained using the Multireference Configu- 
ration Interaction method (MRCI) restricted to single 
and double excitations with a large active space. Scalar 
relativistic effects were included by using the small- 
core fully relativistic energy-consistent pseudopotential 
ECP28MDF [iil from the Stuttgart library. Thus, in 
the present study the Rb2 molecule was treated as a sys- 
tem of effectively 18 electrons. The [lAslApldGf] basis 
set was employed in all calculations. This basis was ob- 
tained by decontracting and augmenting the basis set of 
Ref. [6^ by a set of additional functions improving the 
accuracy of the atomic excitation energies of the rubid- 
ium atom with respect to the NIST database [1^. The 
DEA-IH-FS-CCSD calculations were done with the code 
based on the ACES II program system [70|, while the 
MRCI calculations were performed with the MOLPRO 
code [71]. 

The rovibrational levels of the A^I]+ and b'^n„ ex- 
cited state manifold are computed by diagonalizing the 
Hamiltonian (jl2p represented on a mapped Fourier grid, 
employing about Nji = 512 radial grid points. For the 
calculations in the field we complement our Fourier grid 
representation for the radial part by a basis set expan- 
sion in terms of Legendre polynomials for the angular 
part, taking advantage of the magnetic quantum number 
m being conserved. We find that jmax — 19 is sufficient 
to obtain converged results for I < 2 x 10^ W/cm^. Pres- 
ence of an intense non-resonant field leads to strong hy- 
bridization of the rovibrational levels, and an adiabatic 
separation of rotational and vibrational motion is not 
applicable [H, HI]. We account for this fact by diago- 
nalizing the full two-dimensional Hamiltonian, Eq. (fTOt . 
For I 7^ 0, the non- resonant field mixes different partial 
waves, and j and j' are not good quantum numbers any- 
more. For the sake of simplicity, we label the field-dressed 
rovibrational levels by the field-free quantum numbers. 
Note that the field-dressed levels are adiabatically con- 
nected to their field-free counterparts even for very large 
intensities. 



V. NUMERICAL RESULTS AND DISCUSSION 
A. Potential energy curves 

To test the ability of the ab initio approach adopted 
in the present work to reproduce the experimental data, 
we first check the accuracy of the atomic results. In Ta- 
ble |T] we report the excitation energies at the dissociation 
limit computed with the DEA-IH-FS-CCSD method and 
compare the results to the experimental non-relativistic 
excitation energies obtained from the Lande rule. In- 
spection of Table U shows that the agreement between 
the theoretical and experimental excitation energies is 
perfect. For the 5s -|- ns and 5s -I- np dissociation limits, 
the RMSD is only 21 cm~^, which represents an error of 

0. 08%. When the D states are included this good agree- 
ment is somewhat degraded. The RMSD is now 49 cm~^, 

1. e., 0.26%. This is due to the lack of g symmetry func- 
tions in the basis set used in our calculations. Note par- 
enthetically that we could not include g functions in the 
basis, because the ACESS II program does not support 
g orbitals in the calculations involving pseudopotentials. 
Our method reproduces very well the electron affinity of 
the Rb atom, 3893 cm~^ on the theory side vs. 3919 cm^^ 
measured in Ref. (72j . as well as the ionization potential, 
33630 cm~i vs. 33690 cm"! Finally, we note that 
the ground state static electric dipole polarizability of the 
atom obtained from our molecular calculations is 319.5 
compared to 318.6 ajj from the most sophisticated atomic 
calculations by Derevianko et al. [TSj . 

The computed potential energy curves are reported in 
Fig. [T] for the and •^E+ symmetries, in Fig. [2] for 
the and symmetries, in Figs. [3] and |4] for the 
^Hg and '^Hg, and ^H„ and ■^H„ symmetries, respectively. 
Finally Fig. [S] shows the potential energy curves for the 
singlet and triplet gerade and ungerade states of A sym- 
metry. The spectroscopic characteristics of the singlet 
gerade states are reported in Table iHl while Table |IlT] col- 
lects these properties for the triplet gerade states. Ta- 
bles IIVI and |V] present the results for the singlet and 
triplet states of ungerade symmetry, respectively. Inspec- 
tion of Figs.[T]to[5]reveals that almost all potential energy 
curves show a smooth behavior with well defined minima. 
Some higher states display perturbations, mostly in the 
form of avoided crossings, due to the interaction with 
other electronic states of the same symmetry that are 
located nearby. At high energies the density of states be- 
comes so high that the avoided crossings produce some 
irregularities in the curves. This is especially true for the 
singlet and triplet gerade and ungerade states of sym- 
metry. The H states show less perturbations, except for 
the avoided crossings between the curves corresponding 
to the (3)iHg and (4)ing, and {Sfllg and (4)3Hg states. 
Interestingly, the H„ states and the A states do not show 
any irregularity due to nonadiabatic interactions between 
the states. 

The agreement of the present potentials with those de- 
rived from the experimental data is very good. This is 



TABLE I: Asymptotic enrgies (in cm ^) and molecular states arising from different states of rubiudium atoms 



asymptote 


enerev 
(present) 


enerscv 
(exD.l 


molecular states 


















12731 


12737 


^n,, 

3v+ 3tt 

1 -'--'-11 


^E+ 


(l)'^S(5s) + (l)'^D(4:d) 


19471 


19355 


3v^+ 3^^+ 3-pr 3-rT 
^g I ) -l-lsi -l-'-ii) 


'Ag, 
^g ; 


fl')^5'('5s~)+('2~)^5'f6s) 


20126 


20133 


1^+ lV+ 3^+ 3^+ 
^g 1 5 ^g > -^u 




(l)^S(5s)+(2)2p(6p) 


23732 


23767 


3t+ 3tt 


3V+ 3tt 
^g , iig, 


(l)2p(5p)+(l)2p(5p) 


25462 


25475 


"'"^^(2), ^Tj^ , ^Tlg, ^Hu, ^Ag, 

^^it (2), , "^Hg, ^IIu, 


{lfS{5s)+{2fD{5d) 


25736 


25707 


ly^+ ly+ Itt Itt 

3y^+ 3y^+ St-t St-t 
^g ! ! ^^9: Llu: 


'Ag, 1A„, 
'Ag, 3A„ 




29741 


29771 


^9 1 





TABLE II: Spectroscopic characteristics of the non-relativistic ^|A|g electronic states of *^Rb2 molecule. 



state 


Ref. 


Re 


De 




Te asymptote 






(bohr) (. 


cm-i) 


cm"""^) 






XlEg+ 


present 


7.99 


3912 


56.1 





5s + 5s 




[27] (exp.) 


7.96 


3994 


57.8 











7.96 


3905 


58.4 









present 


10.29 


3102 


32.0 


13545 


5s + 5p 




[36] (exp.) 


10.28 


2963 


31.5 


13602 






m 


10.17 


3084 


31.2 


13559 






present 


10.32 


4210 


32.9 


19180 


5s + 4d 




[41] 


10.20 


4072 


31.9 


19189 






present 


9.34 


4144 


62.0 


19898 


5s + 6s 




present 


9.21 


3483 


37.8 


24166 


5s + 6p 


2nd. min. 


present 


22.22 


2968 


11.0 


24681 






present 


8.93 


3055 


46.6 


24594 


5s + 5p 


2nd. min. 


present 


12.02 


2056 


50.6 


25593 




3rd. min 


present 


34.60 


86 


4.7 


27734 






present 


11.26 


1852 


92.8 


25797 


5s + 5p 




present 


9.47 


183 


41.3 


27465 


5s + 5d 


(i)'ng 


present 


10.25 


1230 


21.7 


15417 


5s + 5p 




[35] (exp.) 


10.24 


1290 


22.3 


15510 








10.24 


1198 


22.0 


15545 




(2)^ng 


present 


9.92 


1326 


31.0 


22063 


5s + Ad 




[41] 


9.88 


1238 


22.0 


22023 




(3)^ng 


present 


9.25 


2833 


43.1 


22149 


5s + 6p 


(4)^ng 


present 


9.48 


2598 


37.1 


22099 


5p + 5p 


(5)^ng 


present 


9.13 


1994 


42.9 


22187 


5s + 5d 


(l)'Ag 


present 


8.18 


5026 


48.7 


18449 


5s + Ad 




[41] 


8.14 


4871 


50.5 


18390 




(2)^Ag 


present 


8.76 


5291 


57.6 


24165 


5p + 5p 


(3)^Ag 


present 


9.22 


2528 


56.5 


27212 


5s + 5d 
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TABLE III: Spectroscopic characteristics of the non-relativistic ^\A\g electronic states of *^Rb2 molecule. 



state 


Ref. 


Re 


De 


We 


Te asymptote 






(bohr) ( 


cm-i) ( 


cm-1) 


(cm-i) 






present 


9.91 


3367 


37.8 


13279 


5s -1- 5p 




m 


9.73 


3345 


36.6 


13298 




(2)^E+ 


present 


8.58 


5372 


51.1 


18017 


5s -|-4d 




m 


8.47 


5347 


51.5 


17914 




(3)^S+ 


present 


9.31 


1657 


38.2 


22384 


5s -1- 6s 


(4)^S+ present 


8.95 


3335 


46.7 


24313 


5s -1- 6p 




present 


9.72 


3488 


19.4 


26065 


5p + 5p 


(6)^S+ 


present 


9.19 


3292 


43.8 


26953 


5s + 5p 


(7)«E+ 


present 


9.12 


3268 


38.5 


27832 


5s -1- 5d 


(i)'n. 


present 


9.54 


-267 


30.3 


16914 


5s -1- 5p 




m 


9.47 


-268 


30.3 


16911 




(2)'n, 


present 


10.56 


3104 


34.2 


20285 


5s -|-4d 




m 


10.53 


2927 


33.6 


20334 




(3)^n, 


present 


9.08 


3416 


45.4 


24232 


5s -1- 6p 


(4)'n, 


present 


9.06 


2646 


27.4 


26735 


5p + 5p 


(5)'n, 


present 


9.09 


2170 


45.8 


27484 


5s + 5d 




present 


8.36 


4181 


48.3 


19284 


5s + 4d 




m 


8.31 


4017 


48.9 


19244 






present 


8.85 


5152 


46.2 


24588 


5s + 5d 




FIG. 1: Potential energy curves for the ^Eg and ^E^ states FIG. 2; Potential energy curves for the ^EJ and ^E+ states 
of the Rb2 molecule. of the Rb2 molecule. 
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TABLE IV: Spectroscopic characteristics of the non-relativistic ^|A|„ electronic states of *'^Rb2 molecule. 



state 


Ref. 


Re 


De 




Te asymptote 






(bohr) 


(cm-i) ( 


cm-^) 


(cm-^) 






present 


9.24 


5967 


44.1 


10680 


5 s + 5p 




m (exp.) 


9.21 


5981 


44.6 


10750 






m 


9.20 


5896 


44.4 


10747 






present 


10.21 


3128 


20.5 


20261 


5s + 4d 




[41] 


10.09 


3003 


22.1 


20258 




2nd. min 


present 


14.11 


3112 


13.5 


20277 






m 


13.81 


2926 


11.5 


20335 






present 


9.37 


1737 


42.4 


22305 


5s + 6s 




present 


9.46 


2390 


31.3 


25258 


5s + 6p 


2nd. min. 


present 


12.64 


2702 


24.3 


24946 




3rd. min 


present 


22.26 


2973 


10.7 


24675 






present 


9.28 


3565 


39.1 


26088 


5p + 5p 


2nd. min 


present 


34.69 


1920 


5.0 


27733 




(6)^2+ 


present 


10.38 


3308 


52.9 


26937 


5s + 5d 


(i)^n„ 


present 


8.57 


1971 


46.9 


14676 


5s + 5p 




[37] (exp.) 


- 


1907 


47.5 


14666 






[41J 


8.48 


1989 


47.9 


14654 




(2)in„ 


present 


8.92 


2369 


31.6 


21021 


5s + 4d 




[37] (exp.) 


- 


2454 


36.4 


20895 








8.77 


2157 


36.1 


21104 




(3)in„ 


present 


9.23 


4927 


40.4 


22721 


5s + 6p 


(4)^n„ 


present 


9.03 


4216 


43.1 


25166 


5p + 5p 


(5)^n„ 


present 


10.06 


3189 


31.4 


26465 


5s + 5d 


(1)'A„ 


present 


9.80 


639 


28.0 


22825 


5s + 4d 




[41] 


9.78 


542 


26.9 


22718 




(2)iA„ 


present 


9.31 


3638 


48.1 


25818 


5p + 5p 


(3)iA„ 


present 


9.40 


2630 


34.2 


27110 


5s + 5d 



demonstrated in Tables [TT] to |Vl where we compare the 
potential characteristics with the available experimental 



data and with the most recent calculations [4l|. For all 
the experimentally observed states, the RMSD of our cal- 
culation is only 75.9 cm~^, i.e., the error is 3.2% on aver- 
age, better than the most recent calculations by AUouche 
and Aubert-Frecon [4l| with a RMSD of 129 cm^^ cor- 
responding to an average error of 5.5%. It is gratifying 
to observe that we reproduce low lying and highly ex- 
cited electronic states equally well. This is in a sharp 
contrast to Ref. [4l| which reproduces the well depth of 
the (2)^n„ state only with an error of 12% compared 
to 3.5% for our calculation. Such a good agreement be- 
tween theory and experiment for the highest observed 
excited electronic state gives us confidence that our pre- 
dictions for the photoassociative production of ultracold 
Rb2 molecules in even higher electronic states 0| are ac- 
curate. Tables HIl to IVl also report the fundamental vibra- 
tional frequencies We for all electronic states considered 
in the present paper. Except for the ground state, the 
agreement between theory and experiment is within a few 



tenths of a wavenumber. Similar agreement was found in 
the calculations by Allouche and Aubert-Frecon [jlj . 



B. Non-adiabatic coupling and spin-orbit coupling 
matrix elements 



The importance of nonadiabatic interactions between 
electronic states, resulting in the avoided crossings of the 
corresponding potential energy curves observed in Figs.[T] 
to O can nicely be explained by analysing the nonadia- 
batic coupling matrix elements computed according to 
Eq. (1181) . The nonadiabatic coupling matrix elements 
are reported in Fig. [5] for singlet and triplet states of E+ 
and symmetry (top) and the 11 states (bottom). As 
expected, the nonadiabatic coupling matrix elements are 
smooth, Lorenzian-type functions, which, in the limit of 
an inifintely close avoided crossing, become a Dirac d- 
function. The height and width of the curve depends on 
the strength of the interaction. The smaller the width 
and the larger the peak, the stronger is the interac- 
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TABLE V: Spectroscopic characteristics of the non-relativistic ^|A|u electronic states of *'^Rb2 molecule. 



state 


Ref. 


Re 


De 


UJe 


Te asymptote 






(bohr) (( 


cm-i) 


(cm-1) 


(cm-i) 




a'^E^ 


present 


11.46 


250 


13.5 


3662 


5s + 5s 




[30] (exp.) 


11.51 


242 


13.5 


- 








11.45 


237 


13.3 


3669 






present 


repulsive 


- 


- 


- 


5s + 5p 


(3)^s: 


present 


11.02 


2761 


40.0 


20628 


5s + 4d 






10.96 


2646 


40.6 


20614 






present 


10.06 


1340 


43.0 


22701 


5s + 6s 




present 


9.18 


2493 


44.7 


25155 


5s + 6p 


(6)^E+ 


present 


9.29 


3235 


40.9 


26147 


5p + 5p 




present 


9.09 


938 


47.2 


28444 


5p + 5p 




present 


7.91 


6969 


57.2 


9677 


5s + 5p 




m (exp.) 


7.81 


7039 


60.1 


9601 






m 


7.88 


7015 


59.7 


9624 




(2)='n„ 


present 


8.73 


3527 


43.5 


19862 


5s + 4d 




m 


8.60 


3497 


43.3 


19764 




(3)'n„ 


present 


9.28 


5117 


40.0 


22531 


5s + 6p 


(4)^n„ 


present 


8.99 


4189 


43.3 


25193 


5p + 5p 


(5)^n„ 


present 


10.04 


3711 


56.5 


25943 


5s + 5d 


(1)'A„ 


present 


9.83 


719 


27.3 


22746 


5s + 4d 




m 


9.86 


619 


25.8 


22641 




(2)'A„ 


present 


9.30 


3695 


40.7 


25761 


5s + 5d 



tion between the electronic states, and the correspond- 
ing potential energy curves are closer to each other at 
the avoided crossing. It is gratifying to observe that the 
maxima on the nonadiabatic couphng matrix elements 
agree well with the locations of the avoided crossing, and 
this despite the fact that two very different methods were 
used in ab initio calculations. Since the potential energy 
curves were shown to be accurate, cf. the discussion in 
Sec. IV Al we are confident that also the nonadiabatic 
coupling matrix elements are essentially correct. 

Rubidium is a heavy atom and the electronic states 
of the Rb2 molecule show strong couplings due to the 
relativistic spin-orbit interaction. Figure [7] reports the 
spin-orbit coupling matrix elements as a function of the 
interatomic separation. The matrix elements are all rep- 
resented by smooth curves approaching the atomic limit 
at large R. The fine splittings of the atomic states are 
very accurately reproduced by our calculations. For the 
first excited P state, the theoretical splitting between the 
1/2 and 3/2 components is 236.2 cm~^ as compared to 
237.6 cm~^ from the experiment. It is also gratifying 
to observe that our ab initio calculations reproduce very 
well the spin-orbit coupling functions obtained from fit- 
ting analytical functions to high-resolution spectroscopic 
data for the A^E+ and b^n„ manifold of states [s^]. This 
gives us confidence that also perturbations in the molecu- 
lar spectra due to the spin-orbit interaction will correctly 
be reproduced from the present ab initio data. 



C. Electric transition dipole moments and electric 
dipole polarizabilities 

A full characterization of the molecular spectra re- 
quires knowledge of the electric transition dipole mo- 
ments. These were calculated according to Eq. (IT7|) and 
are presented in Fig. [8] for transitions from the X^S]+ 
ground state and in Fig. [9] for transitions from the a'^Sj 
lowest triplet state. The strongest transitions from the 
ground singlet state are those to the A^S+ and (l)^n„ 
states, i.e., to states corresponding to the first excited dis- 
sociation limit. All other transition moments are much 
smaller, suggesting that the corresponding line intensi- 
ties in the spectra will be much weaker. The same is 
true for transitions departing from the a'^E^ state. The 
transition moments do not show a strong dependence on 
R, except at small interatomic separations, and smoothly 
tend to their asymptotic atomic value. 

The static electric dipole polarizabilities for the X^I]+ 
electronic ground state, the a'^Ej state, and the rele- 
vant excited A^E+ and b-^IItj states are presented in 
Fig. [51 They show an overall smooth behavior and 
also tend smoothly to their asymptotic atomic values. 
The interaction-induced variation of the polarizability is 
clearly visible while changing the internuclear distance R. 
It is significant for excited states, especially for the A^E+ 
state for which the isotropic part a reaches 8000 Oq, and 
the anisotropic part Aa reaches 6000 Aq. Such large val- 
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FIG. 5: Potential energy curves for the ^Ag/^ and ^Ag/^ 
states of the Rb2 molecule. 



ues of both the interaction-induced variation of isotropic 
and anisotropic polarizabiUties suggest that the influence 
of the non-resonant laser field on the rovibrational dy- 
namics and transitions between the ground state, 
and the A^S+ and b-^IIu states should be significant even 
at relatively weak field intensities. 

Note parenthetically that the transition moments and 
matrix elements of the spin-orbit coupling also change 
when a DC or non-resonant AC field is applied, but the 
changes induced on the rovibrational spectrum are ex- 
pected to be smaller compared to the effects introduced 
within Eq. Therefore, the investigation of the field- 

induced variation of the transition moments and spin- 
orbit couplings is out of the scope of the present paper. 



D. 



Rovibrational spectra in the A^Ej -|- b^IIu 
manifold without a non-resonant field 



We now compare in more detail the ability of our ah 
initio data to reproduce the fine details of high-resolution 
experiments of Ref. In Fig. ITOl' a') . we report the ah 
initio and empirical potentials for the A^E+ and b^IT^ 
states of Rb2 . Inspection of Fig. [TOTa') shows a very good 
agreement. The ah initio calculations reproduce the well 
depth of the A^E+ state within 14 cm~^ on the overall 
depth of 5981 cm~^, i.e., within 0.2%. The agreement 
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FIG. 6: Non-adiabatic radial coupling matrix elements be- 
tween states of E (a) and 11 (b) symmetry. 



for the b^IItj state is slightly less good. The difference in 
the well depths amounts to 70 cm~^ for the well depth of 
7039 cm~^. This represents an error of roughly 1%. Such 
an agreement between theory and experiment should be 
considered as a very good. Also the crossing of the A^E+ 
and b^n„ potential energy curves is perfectly reproduced. 
Our dynamical calculations predict the level w' = 21 to 
be the first rovibrational level corresponding to the A 
state, see the rotational spacings in panel (b) of Fig. (TU] 
This is one quantum higher than predicted by the ex- 
periment [s^l, but the 70 cm^^ disagreement in the well 
depths fully explains this difference. 

Figure [10] also reports the rotational constants for the 
deeply bound rovibrational levels (panel (c)) and levels 
at the threshold (panel (d)). Inspection of Fig. [TUT c) re- 
veals that theory correctly locates all levels that are not 
perturbed by the spin-orbit interaction, and the first per- 
turbed level. The agreement in the rotational constants 
for the rovibrational levels in the middle of the potential 
well is less good, but note the scale on the axis. Over- 



"T 1 T 



2'z*-l'n 




b) 



10 15 20 
R (bohr) 



5 10 15 20 
R (bohr) 



FIG. 7: Spin-orbit coupling matrix elements between states 
of ungerade (a) and gerade (b) symmetries dissociating into 
^S(5s) + ^P(5 p). Black circles and red squares are experimen- 
tal from Ref. 1341. 



all, we reproduce semi-quantitatively the pattern of the 
rovibrational levels in this region of the potentials. Also 
the oscillations of the rotational constants reflecting the 
perturbations due to the spin-orbit coupling between the 
A^Ej and the b"^n„ states are correctly described. This 
is in accordance with the good agreement between the ab 
initio spin-orbit coupling and the data fitted to the exper- 
iment shown in Fig. [T) The agreement of the rotational 
constants for the rovibrational levels near the threshold is 
almost perfect. This is partly due to the fact that in our 
calculations we have used the best long-range coefficients 
from atomic calculations [t^. However, the correct long- 
range coefficient alone would not be sufhcient to obtain 
such a good agreement between theory and experiment. 
In fact, panel (d) of Fig. [TU] shows that theory very pre- 
cisely locates the repulsive walls of the potentials near 
the zero crossing. This is very gratifying for a theoretical 
calculation as this region of the potential energy curve is 
very difficult to describe with ab initio methods. 



E. Perturbation of the spectra by a non-resonant 
field 



Bound rovibrational levels are strongly affected by a 
non-resonant field [1^. We demonstrate in this section 
that not only are the levels shifted in energy and is their 
rotational motion strongly hybridized, but also, for levels 
in the coupled A-'^S+ and b'^n„ excited state manifold, 
the singlet-triplet composition may be changed. Note 
that the non-resonant field mixes different rotational and 
possibly also vibrational states, and in the presence of the 
field, v,j, v',j' are not good quantum numbers anymore. 
However, for simplicity, we do not distinguish between 
the field-free quantum numbers v,j, v',j' and the corre- 
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FIG. 10: Characteristics of the rovibrational levels for the 
— Ot component of the coupled A^Ej and b^IIu mani- 
fold of states in *^Rb2: (a) present and empirical potential 
energy curves [sj], (&) rotational spacings, and j = 1 rota- 
tional constants for strongly bound levels (c) and close to the 
dissociation limit (d). 



FIG. 8: Electric dipole transition moments: (a) between the 
X^E^ ground state and excited states of and ^IIu symme- 
try and (b) between the a^T,^ lowest triplet state and excited 
states of ^E^ and ^Ilg symmetry. 




5 10 15 20 25 5 10 15 20 



R (bohr) R (bohr) 

FIG. 9: Electric dipole polarizabilities for the electronic 
ground state (left) and the first excited state (right). 



spending field-dressed labels v,j, v',j' [ssj . 

Comparing three different intensities, Fig. [TT] il- 
lustrates the effect of the non-resonant field on the 
transition dipole matrix elements for transitions be- 
tween the X-'^S+ ground state and the A^E+ and 
b'^n„ excited states. The transition dipole matrix 
elements are calculated as rovibrational average of 
Eq. (fT7)) for given field-dressed rovibrational levels, i.e., 

/ X^E+\ 

Efe=AiE+,b3n„ ^ X){R)cose\(p.^.^ " V and 

shown for the X^I]+ state ground level in Fig. [TTJa) and 
a vibrationally highly excited level in Fig. [T1T 6). These 
levels could be studied using molecules in a molecular 
beam (a) or produced by photoassociation (5) [l^l- The 
detuning is computed as Ev',j> — Eyj — {E-2p^^^ — E-2g), 
with E2p_^^^, E2g the field-free energies of the atomic lev- 
els. Inspection of Fig. [TT] reveals that the transitions get 
shifted as expected, due to the decrease of all eigenen- 
ergies in the non-resonant field [52, 53]. Moreoever, the 
transition strengths are strongly modified. This mod- 
ification is analyzed in more detail in Fig. [12] for the 
strongest transition from the X-'^E+ state ground level 
indicated by an arrow in Fig. Illf ai. Due to hybridiza- 
tion of the rotational motion, illustrated in Fig.fT^a) and 
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FIG. 11: Transition dipole matrix elements between the 
X^E^ state ground level (a) and a highly excited level (b) and 
rovibrational levels of the A^Ej and b'^n„ manifold for three 
intensities of the non-resonant field in *^Rb2. The binding 
energy of the field-free X^E^ v — 113, j = level is Et — 
8.3cm-\ 



(6) in terms of the projections of the rovibrational wave- 
functions onto the field-free partial waves, the wavefunc- 
tions consist of contributions from several field-free par- 
tial waves between which transitions are allowed. This 
yields a series of rovibrational lines observed in Fig.[T2lc) 
instead of the single line for v — 0,j = to v' — 52, j' — 1 
in the field-free case. For the largest intensity shown in 
Fig. [Til I — 2 ■ 10^ W/cm^, the transition matrix ele- 
ments for the strongest lines are clearly larger than in 
the field-free case. This is rationalized by an alignment 
of the field-dressed levels in the ground and excited elec- 
tronic states, with (cos^ 6) > 0.73 for I = 2 • 10^ W/cm^. 
Correspondingly, the field-dressed wavefunctions are lo- 
calized in the angular regions 9 close to and tt. As a 
consequence, the field-dressed transitions are larger than 
the field-free ones due to the angular dependence of the 
matrix elements on cos 9 [tII . 

Figure [T3] illustrates the effect of a non-resonant field 
on the transition dipole moments for a weakly bound 
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FIG. 12: Partial wave decomposition of the field-dressed rovi- 
brational wavefunctions for the X^E^ state v = 0, j — 
ground level (a) and the v' = 52, j' = 1 level (b) of the 
A^EJ and b^Ilu manifold in *^Rb2. Also shown are the 
electric dipole transition moments between the X^E^ state 
V — 0,j = ground level and the rotational manifold with 
v' = 52 (c). 
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FIG. 13: Transition dipole matrix elements for a highly 
excited rovibrational level («' = 463, -BJ'-=/'^^ = 8.3 cm" ^) of 
the A^E^ and b^Ilu manifold and highly excited X^E^ state 
levels in *^Rb2. 



level in the excited A^E+ and b^n„ state manifold. This 
level is particularly well-suited for the photoassociative 
production of Rb2 molecules [l^, and the analysis of 
Fig. [13] is motivated by a recent proposal for enhanc- 
ing photoassociation rates using a non- resonant field [53j . 
While the calculations of Ref. [H^ were carried out for 
Sr2, a somewhat smaller, albeit still significant enhance- 
ment of the photoassociation rate of about one order 
of magnitude can be expected for Rb2 The non- 

resonant field will affect the spontaneous decay of the 
photoassociated molecules which is governed by the ma- 
trix elements shown in Fig. 1131 The field-free data rep- 
resents a rotationally resolved equivalent of Fig. 3 of 
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FIG. 14: Singlet component of the coupled excited state rovi- 
brational levels vs non-resonant field intensity with v' the 
field-free vibrational quantum number. Data shown for j = 1 
(the behavior for other j is very similar). 



Ref. (23. The binding energy of 8.3 cm- ^ in Fig. [H cor- 
responds to the ground state level v — 113, cf. Fig. 11116). 
A weak non-resonant field splits the two lines originat- 
ing from the j' = 1 level into several ones, similar to 
Fig. [T^ c). The transition strength for j — is almost 
not affected by the weak field. This behavior is similar to 
what has been observed for transitions between weakly 
bound levels of the strontium dimer For a strong 

non-resonant field, the binding energies are shifted and 
the rotational spacings become comparable to the vibra- 
tional one, but the overall behavior is similar, again in ac- 
cordance with Ref. [531] . This implies that a non-resonant 
field may enhance the photoassociation rate without com- 
promising an efhcient stabilization into bound ground 
state levels by spontaneous emission as it was observed 
in Ref. [13. 

Finally, Fig. [T3] analyzes the interplay of the spin- 
orbit coupling and the interaction with the non-resonant 
field for several of the rovibrational levels of the A^E+ 
and b^IItj manifold studied in Figs. [TlTa). [T^ and [T^l 
Surprisingly, the levels from the middle of the well, 
v' = 52, ... , 56, show a remarkable dependence of the 
singlet-triplet decomposition on the non-resonant field in- 
tensity. On the other hand, the singlet-triplet character 
of weakly bound levels of the A'^E^ and b^XI^i manifold, 
shown here for the representative v' = 463, is hardly af- 
fected. This behavior can be understood by inspection of 
the i?-dependence of the polarizability components and 
the spin-orbit coupling matrix elements, cf. Figs. [S] and 
[7l Weakly bound levels have most of their amplitude at 
internuclear separations larger than i? = 20 ao. The spin- 
orbit coupling is strong at large internuclear separations 
and smaller at intermediate separations, while the oppo- 
site is true for the polarizability components. A large 
dependence of the singlet-triplet character on the non- 
resonant field intensity is expected when the interaction 
energy with the field and the spin-orbit coupling become 



comparable. Due to the i?-dependence of the polarizabil- 
ity, for weakly bound levels this requires field intensities 
in excess of 10^°W/cm^. On the other hand, the more 
deeply bound levels, v' — 52, . . . , 56, have their outer 
turning point near R ~ 12 oq where the polarizability is 
large and the spin-orbit coupling is small. Therefore, in- 
tensities of the order of 10^ W/cm^ yield an interaction 
energy with the field that is comparable to the spin-orbit 
coupling. For example, for 1 • 10^ W/cm^, the Stark shift 
of the levels amounts to about 15 cm"'^. Their vibrational 
spacing, of the order of 20cm~^, is also comparable. The 
interaction with the non-resonant field will then affect the 
singlet-triplet character of a rovibrational wavefunction, 
provided the i?-dependence of polarizabilities differs for 
singlet and triplet states. This is indeed the case, cf. 
Fig. O explaining the changes in the singlet-triplet de- 
composition observed in Fig. 1141 



VI. SUMMARY AND CONCLUSIONS 

In the present work we have investigated how the 
spectroscopy of the Rb2 molecule is affected by apply- 
ing a non-resonant field. Our emphasis has been on 
the manifold of the spin-orbit coupled A^I]+ and b^n„ 
excited states. To this end we have derived the elec- 
tronic Hamiltonian describing the interaction of a di- 
atomic linear molecule with a non-resonant field in gen- 
eral and the Hamiltonian describing the nuclear motion 
in a non-resonant field for the manifold of the coupled 
A^E+ and b^n„ excited electronic states in particular. 
We have employed the double electron attachment inter- 
mediate Hamiltonian Fock space coupled cluster method 
restricted to single and double excitations for all elec- 
tronic states of the Rb2 molecule up to the 5s -I- 5d dis- 
sociation limit of about 26.000 cm^^. The agreement be- 
tween the present results and those fitted to high reso- 
lution spectroscopic data is very good, both for the well 
depths and the vibrational frequencies. The accuracy 
of the present results for the potential energy curves is 
much higher than the previous electronic structure cal- 
culations in Refs. [s^ and slightly better than in the 
most recent study by Allouche and Aubert-Frecon [4l| . 

In order to correctly predict the spectroscopic be- 
havior, we have also calculated the electric transition 
dipole moments, non-adiabatic coupling and spin-orbit 
coupling matrix elements, and static dipole polarizabili- 
ties of Rb2 , using the multireference configuration inter- 
action method. To the best of our knowledge, we have 
reported in this paper the very first calculation of the 
irreducible components of the polarizability tensor as a 
function of R for electronically excited states. For the 
spin-orbit coupled manifold of the A^I]+ and b^H^i ex- 
cited states, we have checked the accuracy of the ab initio 
results with the spectroscopic data. Very good agreement 
was found. 

We have investigated the spectroscopy of Rb2 in its 
rovibronic ground state, corresponding to a molecular 
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beam experiment, as well as in highly excited vibrational 
levels of ground and electronically excited states, typical 
for photoassociation experiments at ultracold tempera- 
tures. In both cases, the spectroscopy is significantly 
altered by a non-resonant field. Specifically, fields of 
the order of lO^W/cm^ are found to split a single rovi- 
brational line into several ones and shift the lines by a 
few cm~^. The splitting is due to rotational hybridiza- 
tion, i.e., the field-dressed wavefunctions are made up of 
several field-free partial waves with comparable contri- 
butions. For strong non-resonant fields, of the order of 
10^ W/cm^, alignment leads to an increase of the transi- 
tion strengths compared to the field-free case, due to lo- 
calization of the rotational wavefunctions in regions close 
to 9 — and tt, and the dependence of the transition 
matrix elements on cos 9 [TJ] . We have also investigated 



the effect of a non-resonant field on the transition ma- 
trix elements that govern stabilization by spontaneous 
emission for photoassociated molecules [23 . Similarly to 
strontium molecules [s^, transitions occur to the same 
vibrational levels as in the field-free case. This implies 
that a non-resonant field may be used to enhance the 
photoassociation rate [s^ without deteriorating stabi- 
lization of the photoassociated molecules into bound lev- 
els of the electronic ground state. Somewhat surprisingly, 
we have found a non-resonant field to significantly modify 
the singlet-triplet character of rovibrational levels in the 
A-'^5]+ and b'^n„ excited state manifold for levels in the 
middle of the potential wells, while weakly bound levels 
remain rather unaffected. We have identified two condi- 
tions for a modification of the singlet-triplet character - 
the interaction energy with the field needs to be compa- 
rable to the spin-orbit coupling and the dependence of 
the polarizability tensor components on the interatomic 
separation must differ for singlet and triplet molecules. 



If both conditions are fulfilled, as was found to be the 
case for Rb2 levels of the A^I]+ and b'^n„ manifold with 
vibrational quantum number around 55, the singlet or 
triplet character of a rovibrational wavefunction can be 
controlled by a non-resonant field. 

An interesting perspective for coherent control arises 
when applying a non-resonant field to degenerate excited 
electronic states. We have shown that, for degenerate 
states, a non-resonant field introduces a coupling between 
different states, 0+ and 2u in the present example. In co- 
herent control based on wavepacket motion, such a cou- 
pling between different states can be used to shape the 
effective potential energy curve governing the wavepacket 
dynamics \E^. Using a non-resonant field, for example 
in the far infrared, comes with the advantage of small 
losses even for strong fields. Non-resonant field control 
of photoassociation rates [s^ or wavepacket dynamics 
based on field-induced resonant coupling [59| represents 
a new twist to manipulating molecules with non-resonant 
fields [Zl. 
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